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Abstract 
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1 Introduction 

This paper is a sequel and a generalization of our previous paper [1] . We consider a Herniitean random matrix model 
where the spectrum of the matrix is constrained on a pre-determined union of intervals as well as being subject to a 
(varying) external potential. In this setting the endpoints of the intervals are called hard-edges and — typically — 
the density of eigenvalues near a hard-edge has power-law of the form {x — xq)^^. If — however — the external 
potential confines the asymptotic spectral density away from the hard-edge, the model is — de facto — independent 
of the location of the hard-edge. 

In our previous paper we described the statistics of the first finitely many eigenvalues that start populating a 
newly-forming spectral band and showed that the predictions based on loop equations of [9] were in fact correct 
[16, 4]. Additionally we improved substantially the error estimate, allowing to treat the transitions occurring when 
the population increases by one. We called the point at which the spectral band is about to emerge the spectral 
outpost, picturing the eigenvalues as the first colonies of a large population (living within the main bands). Pushing 
the same analogy we could call the hard-edge a port which has a sea on one side where the colonies cannot live (i.e. 
where the eigenvalues are forbidden to be). The present paper is in the same spirit as [1], but we want to describe 
a similar first-colonization that occurs precisely at a hard-edge of the spectrum (or port). In the process we will 
improve some of the ingredients that appeared in [1]. 

In order to simplify the setup we will simply assume that the model we are describing consists of a model of 
positive definite Hermitean matrices M S 7i+ (here 7i_|_ denotes the cone of positive definite Hermitean matrices), 
thereby putting the hard-edge at the origin. Of course one could as well consider any subset of H determined by the 
requirement that Sp{M) C i-c. the spectrum lies within a pre-determined union of disjoint intervals. Since 

our analysis is local to one hard-edge there is little loss of generality in restricting ourself to the case of positive 
matrices. The reader will recognize that all the considerations can be extended without much effort to the most 
general situation. 

We thus consider the matrix model with (unnormalized) probability measure 

d^ := (det A//)"e^Tr'-^(A^) dAf , a > -1 , M > 0. (1.1) 

The potential V{x) is a (scalar) function with the properties that it is real-analytic at x = (i.e. at the hard-edge), 
it is bounded from below on and grows faster than ln(l + x"^) at infinity. 

Following our idea in [1] we will use a simplified setup which dispenses us with the necessity of a complicated 
double-scaling limit that amounts to a step-wise (infinitesimal) modification of the potential near the hard-edge. 
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In order to explain the setup in more detail, suppose that the (electrostatic) effective potential [17] (p{x) = 
■^V{x) — g{x) + ^i'^ vanishes at the hard-edge x = as Co x'^ with some Co > (the positivity of Co means that 
the main spectral band is confined away from the hard-edge). We then modify the potential (Sect. 3.1) by adding 
a step-like perturbation of the form 

V{x)^V{x) = V{x)^(^2T^J^ + eN{x)^Xj{x) , ^ ■= I ■ (1-2) 

Here XJ is the characteristic function of a small interval J around the hard-edge. The discontinuous character of 
the perturbation makes no trouble whatsoever, since the characteristic function could be replaced by any smooth 
function that is identically 1 on J and identically vanishing in a slightly larger interval and all the analysis would 
be identical. The real parameter >f determines the strength of the perturbation and the constants are crafted for 
later convenience. The function eN^x), which is independent of >r and or order N~'^ , will be specified in due course. 

Due to Dyson's theorem, instead of studying directly the spectral statistics we focus on the associated orthogonal 
polynomials and the corresponding Christoffel-Darboux kernel, in terms of which all correlation functions of the 
eigenvalues can be written [15]. These will be studied with the Deift-Zhou steepest descent method [6] based on 
the formulation of the relevant Riemann-Hilbert problem as in [13]. 

For >f < the orthogonal polynomials do not exhibit any peculiar behavior in the large N limit; for positive 
values of >c new zeroes of the OPs start appearing near the hard edge. The same transitional phenomenon as >c 
crosses a half-integer is observed naturally, here as in [1]: specifically the normalizations we have chosen are such 
that in the asymptotic regime there are K roots near the hard-edge, where K is the integer nearest to >f. Clearly 
transitions must occur a.t x E N + ^. 

At the level of the random-matrix model it will appear that the first K eigenvalues that are growing from the 
hard-edge start populating the new band subject to the statistic of a finite— size K x K matrix model of positive 
matrices subject to the "microscopic" (denoted by the subscript "m" hereafter) measure 

d^i,r^ det(M„0"e-^'- (1.3) 

where Mm is a, K x K positive definite Hermitean matrix (microscopic compared to the N oo original model) 
and VmiC) will be an arbitrary monic polynomial of degree v in ( — a microscopic coordinate. 

Although the main results will be along the same line of [1] there are a few interesting differences. 

• The microscopic potential Vmix) can be directly controlled by e^ix) (1.2). In [1] Vm{x) was only a monomial. 

• The global parametrix requires a new piece called Szego function. It can be written in an arbitrary higher 
genus using Theta functions. 

• The method used in [1] of the "partial Schlesinger transform" was not completely general but seemed to 
require some (upper-triangular) structure of the local parametrix. Here we show that the same method can 
be applied in the most generic situation. 

• The above method provides us with a recursive procedure to construct the asymptotic solution up to an 
arbitrary order of accuracy limited only by the error terms coming from the simple end-points of the spectrum. 

In the main body of the paper we make the simplifying assumption that the support of the equilibrium measure 
consists of one interval (one— cut assumption) confined to the interior of the positive real axis. Of course, this 
assumption is not crucial to the main result. In Appendix C we show (in a somewhat sketchy form) how to 
generalize to an arbitrary number of cuts. 

^The "g-function" is a rather common object in this area of research and will be explained in due course. 
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2 General setting 

2.1 Riemann-Hilbert problem of orthogonal polynomial 

The setting of the present paper is basically identical to [6], [3], and especially to [1]. 
Consider a Hermitean matrix model with measure given by 

^(detM)"e-^*''^(*^MM, 



(2.1) 



where Zn is the normalization constant. The potential V{x) is assumed to be real-analytic on the half-line x > 0. 
In particular V{x) is analytic at a; = 0. 

Let {pn{x)\n — 0, 1, 2, ...} be the corresponding (monic) OPs that satisfy the following orthogonality 



/ pn{x)p„i{x)x°'e T^i'^) dx ^ hnS„ 



The spectral statistics of the model can be computed in terms of the ChristofFel-Darboux kernel [15] 

n-l 



r^i '\ ) Pnix)p„^i{x ) - Pn-l{x)pn{x ) 

K{x,x) ^ 2^ = 

j=0 ^ 



hn-l{x - X') 



(2.2) 



(2.3) 



The asymptotic analysis hinges on the following characterization for the OPs in terms of the Riemann-Hilbert 
Problem (RHP) described hereafter. Define for z S C \ M+ the matrix 



Yiz) := y„(z) 



Pn{z) C[p„](z) 
2-p„_i(z) ^C[pn-l]{z) 



h„-i 



C[p]iz) 



1 f p{x)x°'e T^(^)dx 



2in 



(2.4) 



Here and below we use the coordinate z for the complex plane and x when restricted to the real line. 

The above matrix has the following jump-relations and asymptotic behavior that uniquely characterize it [13, 
14, 11, 12] (we drop the explicit dependence on n for brevity) 



Y+ix)^Y^{x) 



1 2:"e"T^(^) 
1 



on X > , Y{z) - {l + 0{z'^)) 





z^' 



(2.5) 



Replacing the orthogonality condition (2.2) by the above jump (and boundary) conditions (2.5) we obtain the 
Riemann-Hilbert problem for the OPs. Using this setup we especially want to investigate the asymptotics of the 
OPs as their degree n := N + r goes to infinity while r being fixed to an integer. 

2.2 (7-function 

The asymptotic analysis is based on the construction of the so-called g-function; we recall that, in the simple case 
where the contour of integration in (2.2) is the real axis the equilibrium measure is obtained from the solution of 
a variational problem for a functional over probability measures on the real axis, in the sense of potential theory 
[17]. Indeed define the weighted electrostatic energy [17] 



V{x) dfi{x) 



In 



-d/i(a;) d/x(x'), 



(2.6) 



4 




Figure 1: Examples of potentials V and corresponding effective potentials ip for critical hard-edge situations. The 
effective potential behaves as x'^ near the hard-edge, with v — f,2,3,4 respectively (from left to right). The plots 
are numerically accurate, but the vertical axes are scaled differently for V and ip. 

where d/x is a positive measure supported on the positive real axis with total mass T — J^^ dfi{x). 

It is known that the functional attains a unique minimum (under mild assumptions on the growth of V{x) 
at infinity) at a certain measure p{x) that is called the equilibrium measure [17, 7]. 

It is also known [5] that the support of the measure p{x) consists of a finite union of disjoint bounded intervals 
and that p{x) is smooth on the interior of the support. 

Taking avail of the equilibrium measure, the g-function [7] is defined as 

g{z) := J p{x)ln{z- x)dx^Tlnz + 0{z-^), (2.7) 

where the logarithm must be defined with an appropriate cut extending — say — from the leftmost endpoint of the 
support of p{x) to +0O. 

The main properties that enter the steepest descent analysis are the standard properties of the logarithmic 
transform. To this end we note that the representation (2.7) implies immediately that ^g{z) is harmonic away from 
the support of p and continuous on the whole complex plane. The Euler-Lagrange variational equations (equivalent 
to the optimality of the equilibrium measure p [17]) can be rephrased in terms of the following conditions for the 
(j-function. 

• for x e M+ we have 

%(x)>0, ^(z) ^-.g(z) + ^, (2.8) 

for a suitable real constant £. 5R(p is the effective potential of the related electrostatic problem. Since the 
integration (the conductor in the terminology of [17]) terminates at the origin, the variational equations tell 
nothing on the sign of diip on the left of the hard-edge. 

• The opposite inequality (and hence the equahty) holds on the support of p. Especially, £ is chosen such 
that ^(f = on the support of p. (The support of p will be called the cuts because they form the cuts of 
the functions g'{z) and ip'{z).) Here and in the previous point, the g-function should be understood as the 
analytic function defined by its integral representation (2.7) on the simply connected domain obtained by 
removing a half-line starting e.g. at the leftmost endpoint of the support of p and extending towards oo. 
Then ^ip{x) is actually nothing but the boundary- value ^{(p-^-{x) + ip^{x)) on the half- line. 

• In suitable finite left/right neighborhoods of the cuts, the function 5R(^(z), which is also harmonic in the 
domain of analyticity of V{z), is negative. 
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We note that the effective potential is defined within the domain of analyticity of V{z) and in particular in a 
left neighborhood of the hard-edge (the origin) due to our assumption that V{z) is analytic at z = 0. 

In particular ^ip is a real analytic function at a; = as long as the support of the equilibrium measure does not 
contain the hard-edge. If belongs to the support then it is known [5] that — in general — the equilibrium measure 
p{x) has an integrable singularity of type 1 / ^/x. 

We consider the case where ^ supp(p) but the inequality (2.8) is non strict and fails precisely at the hard- 
edge. Figure (2) shows the potential and the effective potential. Taking the cut of g{z) to extend from the 
leftmost endpoints of the support of p to +oo we see that both g{z) and ifiiz) are real-analytic functions in a finite 
neighborhood of the hard-edge x — 0. In particular, because of inequality (2.8) we must have 

ip{x)^Cox''{l + 0{x)) , Co>0, i^eN+, (2.9) 

where v can be any positive integer. The most generic situation is = 1 when, as we will see, the microscopic 
ensemble reduces to a Laguerre ensemble. 

2.3 Local coordinate around the hard-edge 

At the hard-edge the effective potential behaves as ip{z) — ^V{z) — g{z) + ^£ ~ Cqz" with Co > (2.9). We define 
a new conformal parameter z as follows: 

z := Co^'fizy = z + 0{z^) =: ze^'-'K (2.10) 

We define D to be a finite open neighborhood around z — that maps univalently to a disk centered at z = 
by the above relation and such that D n M.^ — [0, e) does not intersect the support of the equilibrium measure. We 
also define the scaling conformal parameter by 

C:=(^— ^(z)j ={CoNr~z where Co (2.11) 

The neighborhood D maps to a disk centered at the origin in the C-plane whose diameter grows as N^. 
We also define the function ri{z), which is holomorphic at the origin, by 



(Co^) or ?=ze''(^). (2.12) 



As can be seen ri(z) measures the mismatch (up to scaling) of the scaling conformal parameter C, with respect to 
the original coordinate z. 

3 Modifying the setting 
3.1 Perturbation of the potential 
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Define J = [0, e) = D n M+; consider tlien the following modified or- 
thogonality relation: 



(3.1) 



where we have defined the exponent of nonregularity 7 := ^ and the 
perturbed potential 



Vix) 



V{x) 
T 



(3.2) 



+ - ( /(C) - (2x + a)7lniV - a7lnCo - a-qix) ]xj{x) 




Figure 2: The potential V with the 
The function /(C) is an arbitrarily chosen real polynomial of degree v-l. chemical potential added, and the corre- 
The combersome form of the perturbation is crafted so as to have sponding effective potential. In this ex- 
the simplest local parametrix and so that the parameter k ^ M. will ample z/ = 1 and the fine-tunmg is thus 
eventually determine the size of the population of the colony near the 
hard-edge. 

The polynomial /(C) is a fine-tuning which allows quite some addi- 
tional generality: this term did not appear in [1] but one could re-read loc.cit. with this fine-tuning in place and 
the obvious changes while maintaining all the results valid. 

Remark 3.1 For the cases v >2 the fine-tuning f{C) is in fact quite a strong perturbation since it scales as N^~'^ 
which is much stronger than the logarithmic perturbation. 



3.2 Normalized and lens-opened RHP 



1 

1 



x = 




Figure 3: The jump matrices for Y. 



In order to streamline the 
derivation we follow [1] and open 
the lenses before normalizing 
the problem, thus modifying the 
jumps as shown in figure 3. We 
thus redefine 



Y 



on the upper lip, 



(3.3) 



Y 



on the lower lip. 



(3.4) 



For the time being the jumps 
on the green circles are the iden- 
tity but later we define separate RHP problems inside the circles. Then we will call all the RHP inside one of the 
green disks the local problem whereas we call the problem outside the outer problem. 
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After the lens-opening we define 



Y{z) 



Ni „ 



N / \ Nf 

r(z)e^T9(^)'^3g2T' 



which satisfies a new RHP: 



e T 



(9- -9+) T-a 



X e T 



gT-(9+-9-) 



X.7(a;) 



Y{z) ~ + 



00 



(3.5) 

(3.6) 
(3.7) 



For simpHcity we assume that there is only one finite band in the spectrum. This means that supp(p) — [a, 6] 
and that the the spectral curve w'^ = (z — a)(z — 6) is of genus 0; the generalization to more bands is not conceptually 
a problem but requires the use of 0-functions which would make the note quite more technical and long. Under 
this assumption, for a; G M, 



g+{x) = —g^{x) + V{x) + ^ , for a; e on the cut, 
g+{x) = g-{x) — 2inT , for x E M_|_ on the right of the cut, 
g+{x) = g-{x) , for X e M+ on the left of the cut. 



(3.8) 
(3.9) 
(3.10) 



Everywhere else on the complex plane g{z) is holomorphic. On account of these properties for the g-function the 
jumps for Y are shown in the figure 4. 

In the following the size of 
the green circles will be fixed to 
a nonzero value. In this case the 
reader can verify that — outside 
of the green disks — the jumps on 
the black and red lines become 
exponentially close to the iden- 
tity, and uniformly so in f] 
L°°. 




x = 



Figure 4: The jump matrices for Y . 



3.3 Outer parametrix: 
Part I 




O 







o 



The variational equations (2.8 
and following) imply that all the 
jumps on the red contours (the 
lips of the lenses) and on the 
spectral gaps are exponentially 
close to the identity jump; after 
removing them we are left with 

the jump matrices as shown in Figure 5. This provides the asymptotic RHP whose solution will be called the outer 
parametrix. 



Figure 5: Jump matrix for k- 
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The outer parametrix ^ will be written in the form 

^^{z):^F{z)^k{z), (3.11) 

where F(z) =1 + Fk,s/z is a unimodular (i.e. det_F'(z) = 1) meromorphic function with a pole only at 2: (at 
the hard-edge); we will call this function a partial Schlesinger transform as in [1]. Remember that the parameter 
>£■ appears in the potential (3.2) and we define the integer K as the nearest integer to >f. Then we define a real 
number 

S -.^ x~ K such that \6\ < ^. (3.12) 

The constant matrix F^.s will be obtained in section 3.6 after we introduce the local parametrix. 
^'if satisfies the following RHP. 

*k(2) ^ (l + 0(z-l))z'■"^ z^oo, (3.13) 



ra-3 

1 


Z 













on the cut. (3.14) 



^K{z) = 0{{z-ayi), ^Kiz) = 0{iz~b)--^), (3.15) 

where a and b are the two turning points. 

In fact, these conditions are still met with the left multiplication by F{z), and we need an additional condition 
to define "ifK- which is the growth condition at the hard-edge. 

*k(z)z~-^'"^ = 0(1). (3.16) 

The above four conditions (3.13)(3.14)(3.15)(3.16) give the Riemann-Hilbert problem for 'i>K- 

3.4 Outer parametrix: Part II (Szego function) 

Let t be the uniformizing map of the genus-0 Riemann surface. We let t^ on the t-plane to map to the hard-edge 
on z-plane. For the simplicity of the normalization we choose the location of the cut and the hard-edge in the 
following way: 

b — a f 1\ b -\- a a — b , ^/l \ 
<t) := ^ + ^) + ^ - - ^0) (- - to) , (3.17) 

where a < b are the endpoints of the band; in the t-plane they correspond to < = ±1. There are 2 points in the 
t-plane projecting to z = namely the two solutions of z{t) = 0. We denote the one outside the unit circle by to, 
the other being 

° to 
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In order to solve a RHP with a;-dependent jumps 
(3.14) we find a function D{x) (usually known as Szego- 
function [18]) that satisfies 



D^{x)D^{x) = x", for X on the cut(s), 
lim D{z) — Doo where Doc G 

z — *oo 

D(x) is analytic and non-zero otherwise. 
It is immediate to find that 



(3.18) 
(3.19) 
(3.20) 



D{x) = 



{t - 1/to) 
t{t - to) 

«/2 



a/2 



i/2 _ 



t - 1/to 



A J y t ' ' (^-^i) 

where the cuts of the non-integer power go entirely in 
the "unphysical sheet" namely between t = and t = 
1/to. Also one can see that -Doo = (^^)"^^- In case of 
several cuts one needs to employ Theta functions and 
we will give the details in an appendix A. 

We now describe a specific solution ^f/f to the above 
RHP (3.13)(3.14)(3.15)(3.16). 

We can write the following solution (a minor modi- 
fication from [1] due to the presence of the Szego func- 
tion) 





I 


\ 










' 


• 

-1 

k 




— • 

1 

to 


1 












Figure 6: The uniformization of the plane sliced along the 
support of the equilibrium measure. The shaded region is 
the "unphysical sheet" . 



^K{z{t)) (-Doo)"^ t'o 



b—a 



l_ 



t-to 

tot-1 



KlT3 



D- 



(3.22) 



Above, t lies outside the unit circle for z not on the cut (the physical sheet), using (3.17). We take the branch 
of the square root such that 



4z'(t) 
b—a 



1 near t = oo. 



3.5 The local parametrix at the hard-edge 

Within the chosen neighborhood of the hard-edge (inside the green circle) the RHP for Y (3.5) reads 

1 2" iv(2>'+")'^e"''(^)e-T(^(^)-2s(^)+£)-/(C) 

1 

1 (^'^N'^'^'fe-^^iQ 



1 ^"7V2><7e-2f v(z)-/(C) 
1 







where K.(C) -^C + fiO- 



We demand that the local parametrix solve the exact same jump within 
i?ji(C), let us first define Rk{C) that solves the following RHP: 

1 (^"e~^'"(^^iV^^^ 







Rk{C)+ = Rk{C) 

Rk{c) - r^"^'o(i) 

Rk{0 -i + o{n-') 



1 



at C G 



when C — !■ 0, 
for some e > on z e 9D. 



(3.23) 

To define the local parametrix 

R+, (3.24) 

(3.25) 
(3.26) 
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Here e is some positive number that will be determined in the subsequent analysis. Increasing e leads to a better 
asymptotics. The second condition (3.25) is equivalent to demanding ^fj^R^ to be analytic, which implies that 
deti? — 1. Note that the singular behavior at the origin for the local parametrix is on the left. 

Note the appearance of the microscopic potential Vm'- in [1] this was completely determined by the effective 
potential and corresponds in the present case to setting the fine-tuning /(C) to zero. However we realize — a 
posteriori — that none of the analyisis in [1] would be affected had we added the fine-tuning, as we are doing here. 

3.5.1 Case >c < ^ 

li >t<\ then the solution to the RHP (3.24) (3.25) (3.26) is written as 



Rk{0 



1 



2111 



C"e 



for X < 



1 

For >f < 0, we may define the local parametrix 

R^iC) := Rk{0 and Fk,s = for x < 0. 

On dD we have 



(3.27) 

(3.28) 
(3.29) 



As we will see, this becomes the error bound of our asymptotics. 

When >c > the error term above will be larger than 0{N~''). We can lower the error bound using the partial 
Schlesinger transform as follows. 

1 ISl^^'f 



C 2iTT 







RK=oiO 



0{N' 



2>t7-27 



for Q < X < 



1 



(3.30) 



3.5.2 Case x>\ 



We recognize in (3.24) the Riemann-Hilbert problem of the orthogonal polynomials for the weig ht e-^'"(«)dC. 
Specifically, if we denote by Pi{C,) the monic orthogonal polynomials that satisfy^ 



dC = m^t 



?7£ > , 



then the solution of the RHP (3.24) (3.25) (3.26) is simply given by 



Rk{0 = z-^'^-^i?K(C), 



(3.31) 



(3.32) 



where 



-2111 



1 

2iTT 

-1 



-Pk-i{0 
Vk-1 Vk-i 



^A--i(e)ec''^"^^^ dg 



^ In the simplest case u = 1 then /(C) is at most an irrelevant constant and the orthogonal polynomials used in the construction of 
the local parametrix are the Laguerrc polynomials with parameter a. 
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^7> 



O 



^K<T3 ^ — 7x0-3 



(3.33) 



It is crucial to point out here that the right multiplier N~''^°'^ is needed to satisfy the correct jump relations (3.24), 
while the left multiplier Cq'^'^'^^ N^^^~^'^'^^ is needed to restore the boundary condition (3.25). 

One can satisfy the boundary condition (3.26) only by choosing K as the closest integer of x\ hence 5 = x — K G 
5) ^ defined already in (3.12). 

We obtain the following estimate holding uniformly on the boundary. 

Rk ^l + 0{N-^+^^^^''), zedD. (3.34) 

We see that if >f e 5 +Z then the error term in (3.34) does not tend to zero (it is 0(1)). This is understandable 
as these values separate regimes where the value of K jumps by one unit and the whole strong asymptotic must 
changes its form. A similar problem arose in [9]. It was shown in [1] and [4] how to overcome the obstacle; there 
is nothing here that differs substantially from that case and the relevant analysis will be recalled later for the sake 
of completeness. 



3.5.3 Improved local parametrix 

To define the local parametrix R^ in this case (and to find the appropriate F(z) appearing in (3.11)), we first 
expand Rk for large C 

/ 



V 



aK 



2iTT a 



'2i7rc27if 



VK-l ^27* 



O 



\ 



where aK is defined from P^iC) — + o-kC," ^ + ■ ■ ■ ■ 
Now we can define the local parametrix R^ as follows. 

Mk.& 



c 



where the matrix Mx,s is a nilpotent matrix defined as follows 



Hk{0 



for X > 



M, 



2ill , Cn 

Vk 



aK 



-aK 



when < 6 < \, 



aK 
2»^ Cg^^ 





VK-l 2 

2z7r ^ 



r2K 



-aK 



when - 5 < (5 < 0, 



when 5 = Q. 



These have been chosen such that R^ gets as close as possible to the identity on the boundary 
that 



R^ = l + 



1 



on 



(3.35) 



(3.36) 



(3.37) 



jY7 mm(l+2|5|,2-2|5|) 

Inside the minimum function, 1 + 2\5\ comes from the term and the 2 — 2\5\ comes from the C"^ term 



We can verify 
(3.38) 
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3.6 Outer parametrix: Part III 

As the last ingredient to the strong asymptotics, let us calculate F{z) := 1 + — |^ (3-11) using the local parametrix 
obtained in the previous section. F{z) is uniquely determined by the following analyticity condition at the outpost. 

0{l) = F{z)^KZ-^'^^(l-^) = (l+^)^KZ-'"'^ (l^^^^—]. (3.39) 

This condition guarantees the analyticity of ^'>ci?jt at the outpost. 
The above condition is equivalent to the two equations below: 

term : = Fk^s [A^ ,Bf] Mk,5 , (3.40) 
z term. U ^ I'kj [A^, ,Ok \^ l^K ■ \ . " ^k,s [A^^ ^^k \ ,p^ , ■ (3-41) 

where the (matrix) coefficients [A^^B^-*] are defined by the expansion below. 

( 1 - =: {aP,bP] + [A(^),bW].+ [a(?),b(^)]z^ + . . . . (3.42) 



Each matrix [A^^B^-*] for j — 0,1, ■ ■ ■ , depends on A'' and has a (finite) limit at — > oo. 

In fact, the equation (3.40) is contained in the other one (3.41). Solving the latter equation for Fk,s and using 
(3.22), we are ready to write down the outer parametrix 

^^{z):= F{z)^k{z) for x>^, (3.43) 

with 

F(z)^l+^ where = [A^^), B^^)] ^ (^[A^^), B^^)] MAi^\ B «] ^) \ (3.44) 
From the behavior of Mk.s (3.37), F{z) is a bounded function of N. 

3.7 Asymptotic solution for Y and error analysis 

Collecting all the results we present the following asymptotic solution for Y (3.5). 

( outside D, 

Yoo ■■= I (3.45) 
[ R,, inside D. 

Though we also need separate parametrices around the edges of the cut, we will skip those because they are well 
known in the literature; for instance, [6] (or our previous [1]). 
To find the error term we define the error matrix as follows. 

£{z) := YooY~^ (3.46) 

The error matrix solves the residual RHP with the jump matrices as shown in the figure 7. It follows from the 
construction that there is no jump inside the green disks, and on the cut. 
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On the disks around the 
edges of the cut (with airy 
parametrix) , the jumps con- 
verge to the identity with uni- 
form error bounds in LP' fl L°° 
as (we simply cite the result, for 
instance from [6]) 

£+£Z^ = 1 + 0{N-^). (3.47) 



The only non-trivial part is 
the jump on which is given 

by 



1 + 0(N-^ mm(l+2|5|,2-2|5|) 




1 + 0{N- 

Figure 7: The jumps of £. 



^7 min(l+2|5|,2-2|5|) / (3.48) 

This is directly obtained from the error bound of (3.38) because i^{z) is uniformly bounded on 9D for growing 
N (3.44). 

A well-known theorem [7] guarantees that error matrix itself is bounded by the same error bound, i.e. £ = 
1 -I- 0{N~^~'^''^^^). This gives the following error term for the (strong) asymptotics of Y . 



^7 inin(l+2|<5|,2-2|5|) 



(3.49) 



4 Roots and kernel at the hard-edge 

4.1 Case |x — K\ ^ \ 

Using the strong asymptotics Y^o (3.45) inside the disk D, let us look at the behavior of the roots at the hard-edge 
in detail. 

The orthogonal polynomial appears at the (ll)-componcnt of the matrix Y . Due to the nature of the transfor- 
mation: Y -^Y (3.5) it is enough to consider Y to study the roots of that polynomial. 

Collecting the previous results, let us write the full expression of Y^o to study the asymptotic behavior for large 

N. 



Y = R 



K.S 



Pk{C) 
-2iTi , , 



(4.1) 



The asterisks denote expressions which are irrelevant to the present analysis. 

Defining a short-hand notation Mk^s ■= (CoNy^MK^s, which is of the order 0(iV2|5|7-7)^ 

we may evaluate 

the leading terms for large N as 



Y — 



C independent 



Pk{C) 

Pk -liC) * 



(4.2) 
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The expression above implies that the microscopic zeroes are asymptotically close to those of Pk{C) li 5 > Q, or 
Pk-i{C,) if (5 < with one stray zero elsewhere. The details on the subleading behavior of the zeroes are basically 
equivalent to the description in [1], section 4.4. Here we summarize the scaling behavior of zeroes when x is not a 
half integer. 

• K < x<K+\: kWK zeroes are at a distance of order 0{N '^'^^) (in the C-coordinate) from those of Pk{C)- 

• K — ^ < X < K . There are K zeroes in total; K — 1 zeroes are within a 0{N~^'^^^^) distance from those of 
of Pk-i{C) Bind the other zero is rimning away to the infinity at the rate ^ N^fl^l in the C-coordinate. We 
may locate this zero as 



Ca 



At 



iV^Tl^l ^here [A'^k^'^P] 



,(0) 



B 



(0) 



(4.3) 



Using the explicit representation (3.22) we obtain 

Considering the relative location a < b and to < (see Figure 6), the above quantity is positive, which 
means that the stray zero runs away to C = +oo (namely towards the main spectral band). Note that in the 
x-coordinate the stray zero actually moves towards the hard-edge, but at a slower rate than the other ones, 
namely as ©(AT-'^^l-'l-i)). 



• X — K: There are K zeroes are at a distance of order 0{N^'^) from the zeroes of Pk{C) ° ^^„] Pk-i{C)- 



The Christoffel-Darboux kernel (in terms of which all the spectral statistics can be computed) is given by 

det 



1 p„(z)p„_i(z') -p„_i(z)p„(z') 



hn-1 



z — z' 

1 



2i7r 



■det 



1 

2i7r" 
Yy^{z) rn(z') 

>2l(2) ^21 (2') 



z' ~ z 



(4.5) 



From the above and from (4.2) we can evaluate the leading behavior of the kernel in the scaling coordinate C- Using 
our asymptotics Y ~ y^Px, we evaluate the kernel in the local coordinate near the hard-edge as 



This implies that the statistics of eigenvalues near the hard-edge in the scaling parameter ^ is governed by the 
microscopic (finite-size) matrix model of "Laguerre" type. 
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4.2 Case \k - K\ 



When K is a, half-integer, i.e. when \S\ — ^, the ^-independent error term (underbraced part) of (4.2) becomes 
0(1). The following fact is useful: 



Mk,s 



= < 



2iTT 
VK-l 



-1] 



-a 



e)(iv-T) 

0(iV"T) 





■ 


1 " 


= -UK 













where we have defined (similar to (4.10) in [1]) 



UK 



2in 





1 



0{N-'<) when(5=i, 
+ O(iV-T) when S = - 



(4.7) 



-a 



j{2K-l) 




By explicitly writing the components as 
the leading behavior of ±1 (3.44) is written as 



B 



(i) 



(i) 



if, 



1 + UK det[Af,A'-^'^] 
iK-idet[B^^\BP] 







(0) 



1 



B 



(0) 

1 









" A (0) 


-a(°) 


a(o) 
2 




A(°) 


-a(°) 







"Bf 


-bI") 


^2 




Bf 


-b(°) 



©(TV-^), 

-o(iv-^). 



(4.8) 
(4.9) 

(4.10) 
(4.11) 



where the subscripts stand for the entries of the corresponding matrix. 
A straightforward calculation gives the first column of Yoo as 



(i^oo)ll 

(^00)21 



Pk{0 



l + UKdet[A^\AW] 



Al' 



(0) 



,(0) 



(4.12) 



-I'-'O 



Pk-i{C) 







\ ^^'^ 1 


b(o) 

2 


+ UK 


. aW _ 



UK det[B^^\AP] 



dct[A(?\A(^)l 



1 + iiK det[A^\A^'] 
]+0{N-^'') for >f=X 



[ aW 1 






. aW _ 


+ 


a(o) 
2 



^A-det[Ag\A^^] 
l + 7.;^det[A(^\A^)] 



a(0) 
2 

1 

2' 



(ijx))ll 

(5^00)21 



£k-iCJA^tPk-i(C) 

(1) T,(0)l 



l+£K-ldet[B^^Bi^^] [ 



+^A-(C) 



b( 



b(°) 
b(°) 

(1) 



(4.13) 



det[A 



^,B(^V4-idet[B(^\B(:!)] 



l + £K-idet[B^',B^'] 



(0) 

1 

(0) 



-^A-iC^A-i(C) 





" bW ■ 




" b(°) ■ 


( 


b(i) 

2 




b(o) 
2 



£g_idet[Bg\Bg)] 
l + i?K-idet[B^\B^)], 



-27N 



for X 



2 



Here, both approximations consist of three terms which — we must note — are each accurate up to multiplicative 
errors of 1 + 0{N^"'). A term containing other dependency on ( — such as ('^Pk{C), for instance — comes with 
0{N~^''') as noted in the formula. 
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The above expressions (4.12)(4.13) can tell the (asymptotic) location of K (or if — 1 for the latter formula) 
zeroes. However, these expressions are not enough to determine the stray zero; for this, we need all the subleading 
terms. Or, the location of stray zeroes have to be computed using the exact improved outer parametrix. 

Finally, using (4. 12) (4. 13) we write the kernel a.t >c ~ K ± ^. 

KniCX ) ^ m-i{C~C') + ^^^(^)^^(^ V ^ ^ 

UKdet[A^\A^^] ^ 
aK ■■= ^ for x = K+\, 4.15 

C^7efte(.)+s(.')-^) ( Pk{Q)Pk-i{C') ~ Pk{(:')Pk-i{0 Pk-i 



K.iCa « r .i:(C-CO - ^^^P^-^iOP.-.iO ] (4.16) 



Pk^i :^ ' ^ 0, for ^ = K-l (4.17) 

l + £K_idet[B«,B(^)] 



These can be considered as "transitional" kernels when the size of the microscopic matrix model is a half-integer. 
Though these kernels do not seem to come from any simple microscopic matrix model, their structures suggest 
very natural interpretation: a mixture (sum) of two adjacent kernels. Such behavior of the kernel has been already 
observed in [1]. 

We conclude this section with the remark that the above analysis requires non-vanishing denominators 1 -|- 
Uif det[A^'', A^'] and l-f detp^'', B^-*]. Using the explicit representation (3.22) we check (below) that the 
denominator does not vanish at least for genus 0. 

.^det[A<^\AW]= _^ f'J" ^ 1-i) jj—^h^) >0 x = X+i,(4.18) 



._.det[B«,B(?)] . ':^l-JZT7^r fl - ^V" ( ^,_JL,.S '^ > for .^K- i.(4.19) 



2^C'j(^^+^'(fe-a)" V toj \ib-a){l-tly^J " 

4-l+2K+a^^_^t2r [ 1^ (6 - a) ( 1 - ig) 

A Szego function in arbitrary genus 

We only sketch the construction since the details would require a good deal of notations to be set up. We will use 
the same notations and ideas contained in [2]. 

We denote by w double-cover of the z-plane branched at the endpoints of the support of the equilibrium measure 

25+2 

w'-.^Hiz^a,) (A.l) 

This is a hyperelliptic algebraic curve of genus g. We denote by oo± the two points above z = c» in the usual 
compactification of the curve, and by p± the two points projecting to the location of the hard-edge^ . We denote 
by LOj the first-kind differentials normalized along the a-cycles: explicitly 

^ (iz 

ujj{z) = aje , (A. 2) 

w 



""The point oo+ is characterized by ui > as 2 S M-|_ near oo. The point p+ is the point on the Riemann surface of w obtained by 
analytic continuation of w on the complex plane slit along (a2j — i,02j) 



17 



where the summation over repeated indices is understood (and they range from 1 to g) and aji is an invertible 
matrix such that loj — djk- 

The notation (rather standard) is lifted from [2] and [10]: the Abel map is understood when writing points as 
arguments of Q and it is based at one of the Weierstrass points (for example ai) A is an arbitrary odd non-singular 
half-period. Recall that (pag. 23 of [10]) all such characteristics A are in one-to-one correspondence with partitions 
of the Weierstrass points into g ~ I and 5 + 3 points {ak^ , • ■ ■ , ctkg^i} U {<^ki ' ■ • ■ ' ^kg+3^ 

Moreover, straightforward computations show that (with some overall ambiguity of signs) 

QAp-^±) - T-diejO)aig (A.3) 
The Szego function we want to define now is a generalization of RHP (3. 18) (3. 19) (3.20). Let 

m 

v{z):=^{z-c,f''^ (A.4) 
i=i 

where pj are some real numbers and the points Cj G M do not belong to the spectral bands. This function is defined 
on a simply connected domain obtained by removing (for example) vertical half-lines in the upper half-plane. The 
Szego function is characterized by the following scalar Riemann-Hilbert problem 

-D+(z)£'_(z) = v{z\ for z on the cut(s), (A. 5) 

lim B{z) = Z^oo where Z^oo e K+, (A.6) 

z — >oo 

D(x) is analytic and non-zero otherwise. (A. 7) 

We denote by oo±, ±7j the points on the hyperelliptic curve above z = oo^z — Cj on the two sheets of the covering 
defined in the customary way; we will refer to the sheet of the points oo+j 7j's as the physical sheet. We leave it to 
the reader to verify that the solution can be written as follows 

. n f ^^j^^^li^^^^^)" .(.)V^. (A.8) 
fJl V0a(p - oo+)eA(p - 7j) / 

The first and third conditions (A. 5) (A. 7) are easily checked. The second condition (A.6) can be verified by taking 
the limit of the proposed expression 

\^ di;QjO)creg J eA(oo+ -7j-) 

Note that -since all branchpoints are real- the Theta- function Qa{p^ oo+) is real for any p G M and outside of the 
spectral bands as follows from ([10], Chapter VI, Prop. 6.1). So is the matrix a^i and the derivatives of at the 
origin. Hence the expression above is real and -in case it is negative- we can simply change the overall sign to the 
expression for D. 

The expression above needs to be supplemented by specifying where the cuts are made; for the denominator the 
cuts run on the physical sheet between and oo+, while for the numerator they run similarly but on the second 
sheet. It is simple to verify (due to the fact that 0a (p — q) has a simple zero for p = q) that in fact the analytic 
continuation of D(x) around each cj yields the same germ of holomorphic function and hence D{x) has -in fact- 
no cuts on the physical sheet. 
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B Partial Schlesinger transform and any order asymptotics 



Here we present a procedure to improve the error bound of the asymptotic solution to an arbitrary accuracy. 
AU-order-asymptotics has appeared in [8]. There, the improved asymptotics are modified by a right multiplication 
to the strong asymptotics. This right multiplication can render the jump matrix become as close to the identity 
as one wishes (since the jump also acts to the right) by a certain well-defined recursive procedure, hence enabling 
one to get a solution with any improved error bound. 

Our method is a little difi^erent in that we use a left multiplication on both outer parametrix and local parametrix. 
In the main text, we have already used such method without which the asymptotics would have failed at any half- 
integer >f. Here we explain how this method can be applied to achieve an arbitrarily high order of accuracy in the 
error commited when matching the outer parametrix to the local parametrix at the hard-edge. It should be said 
that there is no practical purpose in improving the error that comes from the hard-edge (or outpost) parametrix to 
better than 0{N^^) since then errors coming in the matching between the local parametrices at the other (simple) 
turning points are already of that order. In principle the same idea could be applied to the endpoints, but the 
details become different and in the interest of conciseness we refrain from such analysis. 

Inside the disk D, the strong asymptotics of Y is the product of the global parametrix and the local parametrix, 
i.e. Yac ~ ^ ^R^. Outside the disk D the strong asymptotics is simply Y^o — ^'^ and, therefore, we have a jump 

on the boundary 9D. 

To have a good asymptotic solution, this jump must be as close to the identity as possible. To this effect 
we will "improve" R^^^ by a suitable left multiplication. This left multiplication will introduce a new singularity. 
Since the strong asymptotics ^ ^R^t must be holomorphic at the outpost, ^ also needs to be modified by a left 
multiplication (that resembles a Schlesinger transformation) so as to cancel the singularity. Below we describe the 
whole procedure. 

Using (3.32) Rk ■= z^^^^Hk = (-'^"'^ (l + 0(1/C)) C^'^^ we can obtain the following expansion. 

= 1+1^ + § + •••, (B.l) 

z Z'^ 

where Yj's are matrices of order 0{N^'^^) from the consideration of z N^^Q. This grading property is crucial 
for the error analysis. 

Since det Hk = 1 both sides of the above equation (B.2) are also uni-modular. In such case, we claim that the 
following factorization can be done 

. (. . f ) (. . f ) (. . t) (. . I) ^ ^ . (X . t ) (x . t) . <B., 

X ^1 4- O ^ ^ foi' arbitrarily large positive integer p, (B.3) 

where MjS (Mjs) are all nilpotent matrices so that each factor is uni-modular. 

The proof proceeds by induction . First , any uni-modular matrix that has an expansion 1 -I- -I- O ( ) must 
have a traceless matrix A. Secondly, any 2 by 2 traceless matrix can be represented as the sum of two nilpotent 

®We will use the same notations R^t for the "improved" local parametrix and ^t^^ for the "improved" global parametrix. 
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matrices, say A — Afp+i + Afp+i. Then the series can be expressed as 



A 



O 



Zp+2 




(B.4) 



This gives the recursive step of increasing p by one in (B.2). Q.E.D. 

With such expansion in mind, we define the improved local parametrix by 



:= z 



_ ~-A'(T3 





zP 



Rk- 



(B.5) 



Here we have used (1 



M 



for any nilpotent M . Using the grading property that we mentioned 
earlier, we obtain R^ = \ + ©(A^-'^fP+i)), which will be the error bound of our improved asymptotics. 

Now we improve the global parametrix such that ^ :^cRx is holomorphic inside D. We recall that k has a 
singularity at the outpost such that is holomorphic at the outpost. We propose that 



= F{z) - Fp{z)Fp{z) ■ ■ ■ F2{z)F2{z) F^{z)F^{z) ^k, 

= l+ -^ + ^H + H — ^ (and similarly for Fj{z)). 

The (recursive) construction follows. To illustrate the method let us write out the asymptotics. 

part 2 



^ R 



(B.6) 
(B.7) 




(B.l 



part 1 

We understand kz~^''^ is holomorphic at the outpost. First we will find Fn that makes "part 1" holomorphic 
at the outpost. Having a holomorphic "part 1", we will then find Fn that makes "part 2" holomorphic. 

Obviously, these steps are recursive and we only need to clarify one generic step: finding j matrices, Fji^ • • • , Fjj, 
that satisfy the following analyticity condition. 



F. 



33 
Z3 



A{z) 1- 



zi 



)-0{z\ 



(B.9) 



where A(^) is a matrix whose entries are analytic at the hard-edge (or outpost) which is determined from the 
previous step in the induction process (e.g. part 1 and part 2 in (B.8)). Denoting the expansion of A as A(z) = 
Aq + Aiz + A2z'^ + • • ■ , and defining a new expansion: 



AI^'1 := A(z) 
we can write the following solution: 
{Fji, • • • , Fjj) — 

( 



Afz + A'i'z 



m j2 . 



(^[i'l,.--,#,).M, 



••• 
••• 

Ao Ai 





^0 



Ax 



Af 



A 



(B.IO) 
(B.ll) 
.(B.12) 



The global parametrix ^'^^ (B.6) is then fully determined by the recursive procedure. 
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C Construction of the outer parametrix for arbitrary number of cuts 



We refer to the corresponding appendix in [1] since the formulae are identical, the only difference being in the 
picture that accompanies the appendix. 
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